Disorder driven destruction of a phase transition in a superconductor 
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We investigate the effects of disorder on a layered superconductor. The clean system is known to 
have a first order phase transition which is clearly identified by a sharp peak in the specific heat. The 
peak is lost abruptly as the strength of the disorder is increased. Hence, for strong disorder there is 
no phase transition as a function of temperature but merely a crossover which is still detectable in 
the IV characteristic. 
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Disorder is known to play an important role in the 
phase diagram of even so-called clean superconductors 
at low temperatures and high magnetic fields. It has re- 
cently been show that these clean systems have a first 
order transition. It is believed that this transition is as- 
sociated with melting or decoupling of the vortex system. 
The thermodynamic signal of the transition is lost at a 
'critical point' below which pinning is thought to domi- 
nate the behaviour 

Significant effort is currently being invested in at- 
tempting to understand theoretically the effect of dis- 
order on the behaviour of the magnetic flux system in 
superconductors [§JJ]. In this letter we discuss simula- 
tions of the phenomenological behaviour and response to 
disorder of a layered vortex system. Our model is de- 
liberately made sufficiently simple that we are able to 
identify the mechanisms behind the effects we observe. 
We believe that our results demonstrate which degrees 
of freedom of the vortex system are necessary in order to 
interpret the real experimental data. 

We present the results of a 3D layered simulation in 
the presence of point disorder. We have included only 
the degrees of freedom associated with vortex lines, an 
approximation we believe to be valid away from the vicin- 
ity of the zero field transition, T c . The vortex lines con- 
sist of stacks of pancake vortices. These stacks are able 
to cross and to decouple when it becomes energetically 
favourable and this is an essential feature of the model. 
The clean system has a first order decoupling transition, 
with an entropy jump comparable to that seen experi- 
mentally away from T c : namely ~ 0.4fcs/pancake/layer. 
However, for strong disorder the first order transition is 
reduced to a gradual crossover. This can be seen in many 
aspects of the behaviour. The pronounced peak in the 
specific heat that is a feature of the clean system is no 
longer present and the diffusion becomes thermally acti- 
vated with an energy scale set by the strength of the inter- 
layer coupling. Thus the decoupling transition is replaced 
by a depinning crossover whose underlying mechanism is 
plastic cutting of the vortex lines. This is confirmed by 
direct measurement of the cutting frequency. Hence a 
vortex glass phase H cannot occur in this system since 



the divergence of the elastic creep barriers at vanishing 
driving force Q will be cut-off by the finite activation 
energy barrier for plastic cutting of the flux lines. We 
analyze both thermodynamic and transport properties. 
The specific heat behaviour demonstrates that above a 
threshold disorder destroys the phase transition of the 
clean system. However, the transport properties are rel- 
atively insensitive to the disorder. This stresses that one 
must be cautious when inferring the existence of phase 
transitions from transport data |7| . The manner in which 
the clean system phase transition is destroyed by the in- 
troduction of disorder has been investigated by Kierfeld 
et al. [|]|| . The destruction of the thermodynamic phase 
transition may in fact take place via a phase transition 
as a function of disorder in our simulation. 

The model is a layered pancake system which demon- 
strates quantitative features similar to those seen in the 
highly anisotropic Bi-22f2 high temperature supercon- 
ductor. The transition is studied as a function of fixed 
magnetic field, B (density of vortices), which is always 
perpendicular to the layers (along the c-axis), and vari- 
able temperature. In order to be able to study the loss 
of the vortex lattice order the vortex position are varied 
continuously as an underlying discrete lattice can lead to 
spurious phases . The disorder is the simplest possible 
— namely, random point pins. 

The system is equilibriated via Langevin dynamics 
with periodic boundary conditions enforced in all direc- 
tions. For simplicity all the temperature dependence of 
the model is introduced via a noise term [|l0| ; additional 
temperature dependencies of the penetration depth and 
other length scales are neglected. We focus purely on 
the vortex lattice aspect of the melting - that is only 
vortex loops representing fluctuations in the positions of 
the flux lines are included. To enable a sufficiently large 
simulation to study 3D effects Gaussian potentials are 
employed for all the interactions. In 2D the true in-plane 
interaction should be the Kq Bessel function but it is 
only reasonable to use this numerically for very short 
penetration depths as otherwise the long range nature of 
the potential results in excessive relaxation times. For a 
layered superconductor the predicted potentials 111) are 
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even longer. That the qualitative behaviour is correctly 
simulated by the Gaussian potentials can be seen by com- 
paring the 2D simulations of Jensen et al. |L2| who used 
the Gaussian potential and Koshelev jl3) who used the 
physical potentials. 

The pancakes in the planes have an in-plane repul- 
sive interaction between them, in our case modelled by 
a Gaussian potential L/" v , = A v exp(— r 2 /£ 2 ), where r is 
the in-plane distance between the vortices, £ v is the in- 
plane vortex range, and A v is the (fixed) strength of the 
vortex potential. 

Across the layers the interaction is more complicated 
as it must include a mechanism for cutting and recon- 
necting the vortices. This is known to be important near 
the melting temperature, in order to allow for the loss of 
long range phase coherence as seen, for example in the 
pseudo-transformer experiments Jlij . It has been shown 
by Clem [[IB) that in order to model the electromagnetic 
interactions across the layers only pair-wise potentials 
are needed. However, Bulaevskii et al. |l6|] have shown 
that including the lowest order terms of the Josephson 
coupling makes three and four-body terms equally nec- 
essary. For weak interlayer coupling we have found that 
the inclusion of the three body term is sufficient. Hence 
the total interaction across the planes is composed of a 
two body attractive interaction (where all r are in-plane) 

Uy l v , (r\ , rjj ) = —Ai exp(— (r\ - r\ ) 2 /£; 2 ) an d a three body 
repulsive interaction. 

rjlW (J J J\ _ 4 qKp -((^-^) 2 +(^-^') 2 +(^'-'-') 2 )/«3 2 b 
u vv<v"[ r ii r ji r k ) — A 3be y ' 

(1) 

Ai and £j are the amplitude and range of the two body 
interaction and similarly A^ and ^3b are the amplitude 
and range of the three body potential. The latter acts by 
excluding three or more pancakes (two in one layer and a 
third in an adjacent layer) from finding their equilibrium 
location to be within a coherence length in the x-y plane. 
The anisotropy is determined by the strength of the inter- 
layer coupling parameter Ai, and is fixed to be 0.2, which 
corresponds to a highly anisotropic system. That is the 
ratio between the tilt and shear moduli is C44/C66 — 0.01 
for layer and vortex spacings pertinent to Bi-2212 in an 
external magnetic field of IT. This is similar to that es- 
timated for Bi-2212 using the elastic modulii as defined 
in Blatter et al. [jy]. The other parameters are: A v = 1 
and A^b = Ai , for the amplitudes of the potentials and 
£ v = 0.6, & = 0.3, = V2& for the ranges in the Gaus- 
sians. Our unit of length is the average spacing between 
the pancakes. 

For this choice of parameters the clean system loses or- 
der in both the a — b plane and the c-direction over a very 
narrow temperature range which is comparable to the ac- 
curacy with which the temperature can be determined in 
the Langevin simulation. There is a pronounced peak in 
the specific heat at this temperature and a Lee-Kostcrlitz 
fl7t binning of the energies yields an associated entropy 



of ~ 0.4fcs /pancake. Furthermore the activation energy 
between the ordered and disordered states grows with 
system size, indicating a first order transition. 

In two dimensions it has been shown by one of the 
authors |l^] that weak disorder with amplitude in the 
range 0.01 < A p < 0.05 has a significant effect on the 
melting temperature (as defined by the onset of diffu- 
sion), Fig. [y. T a b initially decreases (to approximately 
half its clean value) but then increases as the pinning 
energy become dominant, lower inset to Fig. [l|. This is 
due to the weakening by disorder of the shear modulus of 
the system. In three dimensions this effect is still present 
but much weaker (of the order of a percent), see upper 
inset to Fig. [|. At these small values of A p there is no 
difference in the results obtained for a broad range of 
vortex-to-pin ratios: 0.5 < N v /N p < 2. 

Experimentally the transition moves to lower tempera- 
tures as the magnetic field is increased |ij . At these lower 
temperatures less thermal energy is available and hence 
the pinning energies become more significant. To sim- 
ulate this without changing the magnetic field a larger 
value of the pinning potential amplitude, A p = 0.5 is 
used. The results are dramatically different from the 
clean system. Firstly, the onset of the diffusion which in 
the clean system is very sharp becomes a gradual onset, 
and is in fact thermally activated. The activation energy 
is determined by the strength of the interlayer coupling: 
E ~ l.bAi (from investigating Ai — 0.2,0.3 and 0.5 with 
fixed A p = 0.5 and varying numbers of pins). This sug- 
gests that the diffusion is occurring via plastic cutting of 
the vortex lines. The cutting rate can also be measured 
directly. We define the rate to be the the number of times 
a pancake cuts and reconnects to pancakes in adjacent 
layers, in a given time interval. There is a plateau in 
the rate above the pure system decoupling temperature, 
see Fig. || for all pinning strengths. However, the onset 
temperature for cutting decreases as the pinning strength 
increases, in agreement with the diffusion results. 

In the clean system there is a peak in the specific heat, 
which is a thermodynamic indicator of the transition and 
can be compared with the experimental results of Zeldov 
et al. and Schilling et al.. In the presence of weak dis- 
order (A p < 0.1) the peak is essentially unaffected but 
the peak vanishes rapidly with increasing disorder when 
A p ~ 0.2 and drops away into the background, Fig. ||. 
It should be emphasized that when a peak is present, it 
is as sharp as our temperature resolution. The height of 
the peak is difficult to determine due to the strong fluc- 
tuations in the vicinity of the transition. Our results are 
fully consistent with the ideas of Kierfeld et al. Q that 
there may be a phase transition as a function of disorder. 

We have shown that the thermodynamic phase tran- 
sition is destroyed by directly investigating a thermody- 
namic property (the specific heat). We now demonstrate 
that the existence of the transition can be falsely assumed 
by indirect measurements, such as IV characteristics. 

For Ap > 0.2 there is no longer a true phase transition 
but this is not apparent if a diffused distance criteria, say 
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R 2 (t — t s ) — 1 is used to determine melting. From this 
criteria one finds that disorder merely shifts the melting 
temperature and one would misleadingly conclude that 
the melting transition survives the presence of disorder. 
It is, however, clear from the form of the diffusion curves 
that the diffusion mechanism has changed, see Fig. ||. 
Measuring the onset of diffusion is a close analogue of an 
experimental IV characteristic. This is demonstrated in 
Fig. H with log-log plots for both the diffused distance as 
a function of temperature (at fixed time) and the volt- 
age as a function of temperature (at fixed driving force) . 
Within the numerical accuracy the IV characteristic is 
linear as a function of current at small current. The IVs 
correspond to thermally assisted flux flow (TAFF) | pl| . 
By assuming that the induced voltage has an Arrhcnius 
form we find that the activation energy is of the same or- 
der of magnitude as that deduced from the onset of the 
diffusion data. Hence our system does not have a vortex 
glass phase ||. 

To summarize we have presented a simple model which 
demonstrates that the thermodynamic decoupling tran- 
sition is abruptly destroyed by a finite amount of dis- 
order. We believe that this is the explanation for the 
critical point in the Bi-2212 phase diagram of Zeldov et 
at. Furthermore, our study shows that transport data 
can suggest the existence of a transition even when it is 
known that disorder has destroyed the thermodynamic 
phase transition. 
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FIG. 1. Main figure: Diffusion as a function of temperature 
for the two dimensional system with N v = 1024 for different 
pinning strengths, A p . The insets are the 2D and 3D be- 
haviours of the melting temperature (as defined by diffusion) 
as a function of pinning strength. It is clear that the effect is 
much more significant in the 2D case. 



FIG. 2. Number of cutting events/pancake as a function of 
the strength of the disorder. 



FIG. 3. Specific heat as function of temperature for the 
clean (solid line) and with 144 pins of amplitude A p — 0.5 
and range R p — 0.125 (dashed line). Both systems consist of 
144 pancakes in each of 8 layers. 
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FIG. 4. Average of square of moved distance over a fixed 
number of time-steps for different number of pins JVp.The am- 
plitude is A p = 0.5 and range is R p — 0.125 in all cases. The 
system consists of 144 pancakes in each of 8 layers. The in- 
set shows the temperature at which R 2 = 1 for different pin 
densities 



FIG. 5. Plots of the Log(Voltage) and Log(diffused dis- 
tance) versus Log(Temperature) for Ai — 0.2 and A p = 0.5 
indicating an onset temperature for both ~ 0.04. The in- 
sert shows the change in the behaviour of the diffusion curve 
between the clean and strongly pinned systems. 
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